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We formulate an ab initio downfolding scheme for electron-phonon coupled systems. In this 
scheme, we calculate partially renormalized phonon frequencies and electron-phonon coupling, which 
include the screening effects of high-energy electrons, to construct a realistic Hamiltonian consist¬ 
ing of low-energy electron and phonon degrees of freedom. We show that our scheme, which we 
call constrained density-functional perturbation theory (cDFPT), can be implemented by slightly 
modifying the conventional DFPT, which is one of the standard methods to calculate phonon prop¬ 
erties from first principles. Our scheme can be applied to various phonon-related problems, such as 
superconductivity, electron and thermal transport, thermoelectricity, piezoelectricity, dielectricity 
and multiferroicity. We believe that the cDFPT provides a firm basis for the understanding of 
the role of phonons in strongly correlated materials. Here, we apply the scheme to the fullerene 
superconductors and discuss how the realistic low-energy Hamiltonian is constructed. 
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I. INTRODUCTION 

A quantitative description of the strongly correlated 
materials is one of the most challenging goals in con¬ 
densed matter physics. In particular, an accurate treat¬ 
ment for the lattice degrees of freedom in the strongly 
correlated regime is necessary for a description or even a 
prediction of functional materials such as high-transition- 
temperature (high-Tc) superconductors, thermoelectrics, 
piezoelectrics, and multiferroics. However, the interplay 
between strong correlation and electron-phonon coupling 
has yet to be fully understood. For example, the role 
of the electron-phonon interaction in the cuprate su¬ 
perconductors is still controversial.Recently, it has 
been proposed that the electron correlation enhances the 
electron-phonon coupling.The phonon might coop¬ 
erate with plasmons to realize high-Tc superconductiv¬ 
ity. It has been shown that an unusual coopera¬ 
tion between multi-orbital electronic correlation and the 
Jahn-Teller phonons is the essence of high-Tc s-wave su¬ 
perconductivity next to the Mott insulating phase in the 
fullerides. 

In this paper, we propose that a combination of the 
density functional theory (DFT) and model calculations, 
which is one of the most powerful methods to study the 
strongly correlated materials,can also be powerful 
in studying the electron-phonon coupled systems with 
strong electron correlations. This idea relies on the en¬ 
ergy hierarchy in the electronic structure:^® By the strong 
electronic correlation and the electron-phonon coupling, 
the low-energy bands near the Fermi level Ep, which we 
call target bands, may be heavily reconstructed, while the 
structure of the high-energy bands will not change dras¬ 
tically. Furthermore, at a temperature where the low- 
energy phenomena (e.g., superconductivity) emerge, the 


high-energy states are nearly frozen, i.e., they are nearly 
totally occupied or empty. Then, nearly all the exci¬ 
tation processes occur in t-subspace, the subspace which 
the target bands span (for later use, we define r-subspace 
as the rest of the Hilbert space). The most important 
electron-phonon coupling processes are the couplings be¬ 
tween these t-subspace electrons and phonons. There¬ 
fore, the low-energy physical properties are governed by 
the low-energy electrons and the phonons. 

This hierarchical structure allows us to construct the 
following three-stage scheme:^® 

1. Obtain the global energy structure by the DFT and 
define the low-energy subspace. 

2. Trace out the high-energy electron degrees of free¬ 
dom and derive a low-energy effective Hamilto¬ 
nian (downfolding). The degrees of freedom in the 
Hamiltonian consist of the t-subspace electrons and 
the phonons. 

3. Solve the derived model accurately by the model 
calculation method. 

In this scheme, we take into account the material de¬ 
pendence and the high-energy electronic structure by 
the DFT, and the effects of electron correlation and the 
electron-phonon coupling in the low-energy subspace (t- 
subspace) are considered by the model calculation. A key 
step in the scheme is the step 2, i.e., the dowonfolding 
procedure to derive the low-energy Hamiltonian. 

When we restrict ourselves to the electron degrees 
of freedom and forget about the phonons, there has 
been much effort in the development of the downfolding 
scheme. In this case, the low-energy Hamiltonian would 
consist of the electron one-body (hopping) and Coulomb 
interaction terms. By employing a localized basis such as 
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the maximally localized Wannier function,the de¬ 
rived model has a form of the extended Hubbard model. 
The one-body part describes a realistic hopping struc¬ 
ture in the t-subspace. The effective interaction between 
the t-subspace electrons is a partially-screened Coulomb 
interaction. This is because the high-energy electrons, 
which are traced out, gives a renormalization of the 
Coulomb interaction. We refer to it as a “partially”- 
screened interaction because it does not include screen¬ 
ing processes originating from the t-subspace electrons, 
which are not traced out and remain as active degrees of 
freedom. This partial screening is often calculated within 
the constrained random phase approximation (cRPA),^^ 
which considers the screening effect of the high-energy 
electrons within the RPA. 

The downfolding scheme combined with the model- 
calculation method has been successfully applied to e.g., 
iron-based superconductors,cuprates,^°“^^ transi¬ 
tion metal oxides,and organic compounds. 
Based on these successes, many attempts have been 
done to further improve the scheme. For example, there 
have been proposals to improve the one-body part,^®’^® 
and the interaction part."‘®“'^® Nowadays, GW-based (not 
DFT-based) scheme is also intensively studied. 

Despite much effort in the electronic systems, the ab 
initio downfolding scheme for electron-phonon coupled 
systems has not been established. If we include the 
phonon degrees of freedom, the low-energy model ac¬ 
quires the electron-phonon coupling and phonon one- 
body terms in addition to the electron one-body and 
Coulomb interaction terms. As in the case of effec¬ 
tive Coulomb interaction between the t-subspace, the 
electron-phonon coupling and phonon frequencies used in 
the low-energy Hamiltonian should be a partially renor¬ 
malized quantity.®® They are renormalized due to the 
coupling between the phonons and the high-energy elec¬ 
trons. The coupling between the phonons and t-subspace 
electrons is considered when we solve the model by the 
model-calculation method. When we derive the model, 
the renormalization originating from the t-subspace is 
excluded to avoid the double counting of it. Basing 
on this idea, we recently proposed an ah initio scheme, 
called constrained density-functional perturbation theory 
(cDFPT).®! 

In this paper, we elaborate the practical detail of the 
cDFPT method. We show that the cDFPT method can 
be easily implemented by a slight modification of the con¬ 
ventional DFPT method, which is implemented in sev¬ 
eral ab initio packages. Then, we apply the scheme to 
the alkali-doped fullerides,where both the electron 
correlations and the electron-phonon interactions are im¬ 
portant to explain the phase diagram.^® By comparing 
the cDFPT results with the DFPT results, we discuss 
how the partially screened quantities, which are used as 
an input for the model calculation, differ from the fully 
renormalized quantities. 

This paper is organized as follows. In Sec. H, we 
review the DFPT®®^®® to introduce our notation, since 


the cDFPT method is closely related to the conventional 
DFPT. Then, we move onto the the explanation of the 
cDFPT method in Sec. HI. There, we provide practical 
details to implement the cDFPT method. We also briefly 
compare the cDFPT method and another downfolding 
method proposed in Ref. 59. In Sec. IV, we show the 
cDFPT results for the alkali-doped fullerides and com¬ 
pare it with the DFPT results. Finally, in Sec. V, we 
give a summary of the paper. 


II. REVIEW OF DENSITY-FUNCTIONAL 
PERTURBATION THEORY 

The cDFPT method is based on the DFPT method, 
where the fully renormalized electron-phonon coupling 
and phonon frequencies are calculated. As we will show 
below, the cDFPT method can be formulated as a slight 
modification of the DFPT algorithm. Here, we briefly re¬ 
view the DFPT method®®^®® just to introduce a notation 
used in the paper. 


A. Phonon frequencies 

1. Expression for interatomic force constants 


In solids, the phonon frequencies are determined by 
the following equation:®® 

(q) = Wqj.e“(qz/) (1) 

k' Oi' 

with a momentum q, the index for atoms k, and the 
direction of the displacement a = {x,y,z}. This equa¬ 
tion shows that the phonon frequency uiqy is given by 
the square root of the eigenvalues of the dynamical ma¬ 
trix D(q). Since the dimension of the dynamical matrix 
D(q) is 3n with n being the number of the atoms in the 
unit cell, there exist 3n solutions (normal modes), which 
we label by the index v. The eigenvectors of the dynam¬ 
ical matrix satisfy the orthonormality: 




( 2 ) 


The dynamical matrix is related to the interatomic force 
constants (7““/ (q) by 


DZ' (q) = 


y/MZZ. 


■CZ' (q), 


( 3 ) 


where is the mass of the Kth atom. The interatomic 
force constants are written as®® 


cz (q) = 



f ( dp{r) \ 

[1 

\du%{q)) 


-dr 


<- . , d^Vionjr) , 
^ ^ 9n*“(q)an“;(q) 
An 


9u*“(q)a<;(q) 


( 4 ) 
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with the number of the unit cells in the Born-von Kar- 
man boundary condition N, the displacement of the ion 
u, the electron density p, the ionic potential Idon, and the 
Coulomb interaction energy among the nuclei -En- On 
the r.h.s. of Eq. (4), the first (second) term describes 
the contribution from the linear (quadratic) electron- 
phonon coupling and the third term describes the ionic 
contribution. 


2. Electron density response 


with pq being the electron density in the absence of the 
ionic displacement. Eqs. (5) and (7) are the equations to 
determine the electron density response, which are solved 
self-consistently. 

In the DEPT, in order to avoid the cumbersome sum¬ 
mation over the unoccupied states in Eq. (6), one alter¬ 
natively solves the following equations [Eqs. (72) and 
(73) in Ref. 58]: 

('^SCF + Q — £n) I Al/'n) = — (^F,n — E„) AEsCF |' 0 n) 

( 8 ) 


In order to evaluate the interatomic force constants, 
we need to calculate the electron-density response to the 
ionic displacement dp{r)/du^{q}, which is a key quantity 
in the cDEPT method, as we will show below. Before go¬ 
ing into the explanation of the cDFPT, we show how the 
electron-density response is calculated in the usual DEPT 
method. Here, we consider a metallic case.®® In the DFT 
calculation for the metal, it is usual to introduce a smear¬ 
ing function 5(x) and the corresponding smoothed step 
function 9{x) = S(x')dx'. In the present calculation, 

we employ the gaussian smearing S(x) = exp(—x^) 
Then, the expression for the electron density response 
Ap(r) to the ionic displacement is given by 

Ap(r) = ^ - ^'^’"" C(r)V’m(r)(^m|AVscF|^n) 

n,m 

= 2Xl^n(r)A?/’n(r), (5) 

n 

where we define A'(/i„(r) as 
Atpriir) = ^ _ ^^’"" ^m,n'0m(l~)(V^m|APsCF|'0n)- 

m 

( 6 ) 


with composite indices for the band and the momentum 
n,m, the Kohn-Sham (KS) wave function and the 
KS eigenenergy e„. Here, 9p n and are defined as 
0F,n = 0[{eF-Sn)/(T] and 9m,n = ^[(Em-en)/o-], respec¬ 
tively, with the Fermi energy £p and a smearing width 
cr. In the actual calculation, the electron density response 
Ap and the modulation of the potential AVscp have in¬ 
dices of the momentum q, the displaced atom k, and the 
direction a, which we omit for simplicity. The change 
of the potential AVscf due to the ionic displacement is 
given by a sum of the change of the ionic potential Aldon 
and the screening contribution from the Hartree and ex¬ 
change channels (the second and third terms on the r.h.s. 
of the following equation): 


APscF(r)=AMo„(r)+e2 f 
, dVy^cip] 


dp 


Ap(r) 


(7) 


P=Po(r) 


where 

Q — ^ ^ CTm I |; Pn — ^ ^ /5n,m | '^m } | (9) 

m m 

with 


— ^ F.n^n.m H" ^ F.m^m.n Ql- 


Opn. — Of 


’ >771 n 
Ur 




( 10 ) 


Here am’s are parameters to avoid null eigenvalues of the 
(Hscf+Q—£ n) matrix, which can be set to be a constant 
value which is larger than [(maximum energy among par¬ 
tial occupied states) — (minimum energy of the occupied 
states)] for all the partially occupied states, and zero for 
the totally unoccupied states.®® This am parametrization 
enables the calculation without any information about 
the totally unoccupied states. In Appendix A, we show 
that the solution of Eq. (8) is indeed identical with that 
of Eq. (6). 

When the perturbation has the periodicity with the 
lattice (q = 0), the Fermi energy may change and Ap 
acquires an additional term:®® 


Ap(r) = 2 E C(i’)A'!/'„(r)-bp(r,£p)A£F (II) 

n 


with 

= E (12) 

n ^ ' 

The change in the Fermi energy Aep is determined by 
the charge neutrality condition.®® 


B. Electron-phonon coupling 

When the ions move from their equilibrium position, 
the ionic potential changes. Then, the surrounding elec¬ 
trons will respond to the potential change and screen it. 
The electron will feel this screened potential change and 
will be scattered. This process is expressed by the Hamil¬ 
tonian 

7 iel-ph = /— 9n'n0^i 

^ knn'cr 

(13) 
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Here, 


5n'n(k,q) = 




2M,u;q, 

9HscF(r) 


n'k+q 


e“(qi/) X 


^M“(q) 


(14) 


is the electron-phonon-coupliiig matrix element involving 
the Bloch states ijjn'k and ipri'k+q and the ^th branch 
phonon with the wave vector q. annihilates 

(creates) an electron on the nth Bloch orbital with the 
wave vector k and the spin cr. bqi, is the annihilation 
(creation) operator for the phonon labeled by the nth 
branch and the momentum q. 


III. CONSTRAINED DENSITY-FUNCTIONAL 
PERTURBATION THEORY 

A. Basic idea and practical implementation 


Our goal is to derive the low-energy Hamiltonian for 
the electron-phonon coupled systems, which consists of 
the low-energy (t-subspace) electrons and the phonons. 
The Hamiltonian reads 

a = Hel + 'Hel-el + Hel-ph + Hph + ^DCi (15) 

where Hei is the electronic one-body part (onsite energy 
and hopping terms), and "Hei-ei is the Coulomb interac¬ 
tion term, such as the Hubbard U. In this paper, we fo¬ 
cus on the electron-phonon coupling Hei-ph and phonon 
one-body term Hph, which are given by 


the model analysis step. To make it clear that these are 
partially renormalized quantities, we attach the super¬ 
script (p). 

In the following, we discuss how the partially renor¬ 
malized phonon quantities are calculated from first prin¬ 
ciples. For the partial renormalization, we first define 
the bare phonon frequencies and electron-phonon cou¬ 
pling. We then divide the renormalization processes into 
the low-energy contribution, which is to be excluded to 
realize the partial renormalization, and the rest of the 
contribution, which involves the high-energy electrons. 

First, we consider the phonon frequencies. As we see in 
Sec. HA, the interatomic force constants [Eq. (4)], which 
give the phonon frequencies, consist of several contribu¬ 
tions. Since the low-energy Hamiltonian in Eq. (15) has 
the linear electron-phonon coupling term, which gives 
a renormalization of the phonon frequencies, we define 
(ionic contribution) -|- (contribution from the quadratic 
electron-phonon coupling) as “bare” term, and (contribu¬ 
tion from the linear electron-phonon coupling) as “renor¬ 
malizing” term. Then the interatomic force constants 
C““/(q) given in Eq. (4) can be divided as (^^““/(q) = 
^--^^/(q) +^“ Cr/(q), where (q) gives the 

“bare” phonon frequencies 


bare 


CT' (q) = T7 




N[ du%°‘{q)du‘^:{q) 

d^Vionir) 


PiT^) 


^<“(q)5<''(q) 


dr 


,( 18 ) 


and (q) gives the renormalization of the phonon 

frequencies through the linear electron-phonon coupling 


"Hei-ph = ''(k,qXk+qCjk(V + biqu)- 

qL' kz_7 (T 

(16) 

and 

= (17) 

qi/ 

respectively. Here, we employ the Wannier gauge for 
the electronic degrees freedom labeled by i,j, since it is 
convenient for the low-energy solvers to take the Wan¬ 
nier gauge. Hdc is a double counting correction, which 
is discussed in detail in Sec. HID. In this section, we 
show how the phonon frequencies and the electron- 
phonon coupling in the low-energy model should be 
parametrized.®^ As in the case of the effective Coulomb 
interactions in TLei-ei calculated by the cRPA method,^® 
they should be partially renormalized quantities, which 
take into account the renormalization effects associated 
with the elimination of the high-energy degrees of free- 
dom(see Appendix B for the comparison between the 
cDFPT and the cRPA). In other words, we derive the pa¬ 
rameters with avoiding the double counting of the renor¬ 
malization effects which are to be taken into account in 




1 

N 


dp{r) y ^bon(i~) ^^ 

9Ki^)) 5</'(q) 


(19) 


Next, we consider the bare and renormalizing contribu¬ 
tions to electron-phonon coupling [Eq. (14)]. The deriva¬ 
tive of the self-consistent field potential 9 Vscf (r) / du^ (q) 
in Eq. (14) is also decomposed into the bare contribution 


bare 


9FscF(r) 


L 5<(q) J 


5Fion(r) 

a<(q) 


( 20 ) 


and the screening contribution (the change of the Hartree 
and exchange potentials) 


ren. 


'SVscfW 
. '9<(q) 



|r-r'| 

dpjr') 

au“(q) 


dI4c(r) 

dp 

dr'. 


i5(r 


r') 


( 21 ) 


We see that the origin of the renormalization of the 
phonon frequencies and the screening for the electron- 
phonon couplings is the coupling between the lat¬ 
tice and the electrons, and the resulting modulation 
of the electron density due to the lattice displace¬ 
ment 9p(r)/9u“(q). The electron-density modulation 
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dp{v)/du^{q) calculated in the conventional DFPT 
scheme is a sum of the contributions from all the pos¬ 
sible particle-hole excitations [Eq. (5)]. In the cDFPT 
method,®^ we exclude the targetotarget excitation pro¬ 
cesses from the sum in the calculation of the electron- 
density modulation. We use the resulting electron- 
density modulation for the renormalization contributions 
in Eqs. (19) and (21), which are added to the bare con¬ 
tributions in Eqs. (18) and (20). This procedure gives 
the partially renormalized phonon frequencies and the 
electron-phonon couplings. 

Now, we propose a practical way to exclude the 
targetotarget processes from Eqs. (5) and (8), the equa¬ 


tions which determine the change of the electron density. 
If |'0„) in Eq. (8) belongs to the t-subspace, in order 
to exclude the targetotarget polarization processes, the 
r.h.s. of Eq. (8) should be modified as 

(’HsCF + Q — en) I = —Pr{0F,n — E„) AVsCF | '0n) 

( 22 ) 

with Pr being the projection onto the r-subspace. The 
very same constraint can be achieved by solving Eq. (8) 
with modified /3„,m,’s (/3n,m’s) given by 


J 


/3 


n,m — 


dF,n (n, m € t-subspace), 

^F^n^n^m P ^F,m^7n,n P en—Srn ^the Other CaSeS^ . 


(23) 


Note that in the latter case, /3n.m has exactly the same 
form as that of Eq. (10), i.e., Pn,m = /3n,m- Only when 
n,m € t-subspace, l3n,m is modified. We can easily show 
that the r.h.s. of Eq. (22) with the original /3n,m’s is 
equal to that of Eq. (8) with I3n,m’s, which ensures the 
equivalence of the two types of the modifications. Using 
/3n,m in Eq. (23) is also useful to exclude the contribu¬ 
tion to the electron-density modulation from the possible 
change in the Fermi energy in the case of q = 0 [the ad¬ 
ditional contribution given in Eqs. (11) and (12)]. The 
possible change in the Fermi energy originates from the 
intraband transitions at the Fermi level, which are the 
transition processes in the t-subspace and hence are ex¬ 
cluded by employing j3n,m- 

When we consider the practical implementation, if one 
has a code of the conventional DFPT, it is easier to 
modify /3n,m into Pn,m than to employ Eq. (22). One 
has only to modify the part where the I3n,m parameters 
are defined, and no modification is needed in the other 
parts. In Appendix C, we propose an example how we 
modify a source code to introduce I3n,m in the case of 
QUANTUM ESPRESSO package.With Pn,m’s and fol¬ 
lowing the very same flow of calculations of the usual 
DFPT method, one can calculate the electron density re¬ 
sponse to the ionic displacement without targetotarget 
polarization processes. Then, with the resulting electron 
density response, we evaluate the partially-renormalized 
quantities and g^P\ 


B. Relation between fully and partially 
renormalized quantities 

In this section, we show the relation between the 
partially and fully renormalized quantities. The par¬ 


tially (fully) renormalized quantities are calculated by 
the cDFPT (conventional DFPT) method. The electron 
density response Ap to the change of the ionic potential 
AUion (bare perturbation) is given by®^ 

Ap=x°(l-i}x°)-'AVjo„ (24) 

= XDFT 

=X°AUscf, (25) 

where AVscf is the screened potential change, given by 

AUscf = (l-iix°)”'AUon. (26) 

Here, v is given by u = u + K^c with the bare Coulomb 
interaction v and the exchange-correlation kernel iXxc = 
SVxc/Sp (14c is the exchange-correlation potential). Note 
that Eqs. (25) and (26) correspond to Eqs. (5) and (7), 
respectively. The screening expressed in Eq. (26) can 
be divided into two screening steps: One involving the 
high-energy degrees of freedom 

^^SCF = (l-^X°)”'AUon (27) 

and the other associated with the target-target processes 

AUi4 = (l - AuJS^p. (28) 

Here, the total irreducible polarization x° is divided into 
X? and Xr with the polarization within the t-subspace Xt 
and the rest of the polarization Xr = ~ Xt ■ We have 

introduced the superscript p, and / to explicitly distin¬ 
guish between the partially (p) and fully (/) renormal¬ 
ized quantities. is the partially screened Coulomb 

interaction given by 

W<^P) = (l-vx^y^i. 


(29) 
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Since the electron-phonon coupling g represents the scat¬ 
tering of the electrons by AI^cf, the screening process 
for the electron-phonon coupling can be decomposed in 
the very same way as that of AVscf [Eqs. (27) and (28)]; 
that is, = (l — g^^^ is decomposed into 

(30) 

and 

gif) = g(p), (31) 

Eq. (31) tells us that when we take into account the 
target-target screening processes at the DFT level for the 
model with the partially-screened Coulomb and electron- 
phonon interactions, we come back to the fully-screened 
electron-phonon interactions. 

The similar decomposition also applies to the renor¬ 
malization of the phonon frequencies. In this case, 
the phonon self-energy is decomposed. The renormal¬ 
izing contribution to the interatomic force constants in 
Eq. (19) can be recast as 

■■-■C=|g'(')|2xDFT, (32) 

where gi’^) with being the bare 

phonon frequency. For simplicity, we have omitted the 
indices and represent the masses of the nucleus by a sin¬ 
gle mass M. We define the phonon self-energy in the 
DEPT scheme as 

ren. 

The contribution to the phonon self-energy can be di¬ 
vided into Ej and S^, i.e., 

E = Et + E,. (34) 

Here, E^ = \gib)\'^x’DFT with = Xr(l “ «Xr) ^ de¬ 

notes the phonon self-energy due to the electron-phonon 
coupling involving the r-subspace electrons. The other 
part of the self-energy E* = Is'-^^Pxdft with Xdft = 
Xt (l — W^P'^Xt) originates from the coupling between 
the t-subspace electrons and the phonons through the 
partially-screened coupling g^P\ See Appendix D for the 
proof that Et -h E^ = + \gib)\'^^-^^^ is in¬ 
deed identical to E = Xdft • The decomposition 

of E into Et and E^. corresponds to the division of the 
density-response contribution to into the target- 

target contribution and the others, as the cDFPT scheme 
does. With the decomposition of E, we can define the 
partially-dressed phonon Green’s function Di^) as 

[25(p)]-i = [£)(b)]-i _E„ (35) 

with the bare phonon Green’s function Di^). The bare 
phonon frequency is given by the pole of Di^). Simi¬ 
larly, the phonon frequency in the low-energy Hamil¬ 
tonian is given by the pole of DiP). If we further consider 
Ei, we obtain the fully-dressed phonon Green’s function 
Dif) as 

-E*. (36) 


C. Flow of the calculation and practical issues 

As we already mentioned, the flow to the cDFPT calcu¬ 
lation just follows that of the usual DEPT. The difference 
comes from the setting of fin,m parameters. The flow of 
the calculation is as follows: 

1. Optimize the atomic positions within the DFT. 

2. Calculate the global energy structure by the DFT 
for the optimized structure and choose the target 
subspace for which we construct an effective Hamil¬ 
tonian. 

3. Set /?n,m parameters according to Eq. (23). 

4. Perform the phonon calculation with Pn,m parame¬ 
ters (the procedure is the very same as the conven¬ 
tional DEPT case). 

—>■ Obtain the partially renormalized phonon fre¬ 
quencies ujqJ to be used in Eq. (17) and the par- 
tially screened potential change ■ 

5. Take the Wannier matrix element of \ to 
obtain the partially renormalized electron-phonon 
coupling term in Eq. (16) as follows: 



where we use the superscript (w) to make it clear 
that the wavefunction is in the Wannier gauge. 


Finally, we mention one practical issue in obtaining the 
partially screened phonon frequencies . In obtaining 
the fully renormalized phonon frequencies, we often im¬ 
pose the acoustic sum rule to ensure that the frequency 
of the acoustic phonon at q = 0 is zero. To obtain the 
partially renormalized phonon frequencies, we impose the 
same correction of the acoustic sum rule as that used in 
the calculation of the fully renormalized phonon frequen¬ 
cies. Then, the partially renormalized phonon frequency 
of the acoustic phonon at q = 0 does not always go 
to zero. This is because the phonon self-energy involv¬ 
ing t-subspace electrons E* = Is^^^Pxdft can be finite, 
since there can be a finite coupling between the acoustic 
phonon and t-subspace electrons through the umklapp 
(G 0) processes, while the coupling for q-|-G = 0 pro¬ 
cess is zero. We also give another explanation for pos¬ 
sible non-zero phonon frequency for the acoustic mode 
at q = 0 . For example, in the case where the unit cell 
consists of a single atom, at q = 0, the ionic contribu¬ 
tion to the interatomic force constant [the third term on 
the r.h.s. of Eq. (4)] is zero. The first (second) term 
on the r.h.s. of Eq. (4), which is related with the lin¬ 
ear (quadratic) electron-phonon coupling, gives negative 
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(positive) contribution to the interatomic force constant. 
Since the first and second terms cancel with each other, 
the fully-renormalized phonon frequency at q = 0 goes 
to zero. In the cDFPT, we exclude the target contri¬ 
bution to the first term, thus imbalance occurs between 
the hrst and second terms, which makes the partially- 
renormalized phonon frequency nonzero. 

D. Double counting correction 

When we combine the DFT and the model-calculation 
methods, we usually need a double counting correction. 
In the case of our scheme, we have a double counting 
problem for a possible change of the equilibrium positions 
of the atoms due to the coupling between the lattice and 
the t-subspace electrons. The low-energy Hamiltonian 
should be formulated such that we obtain the equilibrium 
positions of the ions which agree with the optimized po¬ 
sitions within the DFT level, after we solve the model at 
the static mean-held (DFT) level. To realize this, we need 
a double counting correction in the low-energy Hamilto¬ 
nian, whose form is 

"^DC = 5lf''(k,q = 0)(c^k^cjk)(&o^ + 

^ kijcr 

(38) 

Here, is the expectation value evaluated within 

the DFT. 

To understand the physical meaning of the double 
counting correction, we consider a simple case, where the 
t-subspace consists of a single band and only one Hol¬ 
stein phonon couples to the electron locally. Then, the 
electron-phonon coupling term in Eq. (16) is given by 

■Hei-ph = (39) 

i 

where we switch to the real space representation and I 
is the site index. ni is the density operator for the site 
I and xi is the displacement of the lattice. The double 
counting correction [Eq. (38)] becomes 

= -^g^^\ni)xi. (40) 

i 

If we put together the electron-phonon coupling, double 
counting, and potential energy terms [the phonon-related 
part of the low-energy Hamiltonian in Eq. (15)], it is 
given by 

I I 

+ ^ i(cc(P))2(a:i - x°)^ + const.(41) 

I I 

where x^ = {ni)/. In the above expression, we 
take the atomic mass to be 1 for simplicity. Now, the 


physical meaning of the double counting correction be¬ 
comes clear: It gives the shift of the potential minimum 
of the lattice vibration by x^, which is proportional to 
the occupation of the electron (ni). x'f gives the equi¬ 
librium position of the lattice vibration without the ef¬ 
fect of the low-energy electron manifold. When we solve 
the model at the mean-field (DET) level, the equilibrium 
position goes back to the optimized position within the 
DFT, since the contributions from Eqs. (39) and (40) 
cancel with each other. 


E. Comparison between our scheme and the 
scheme proposed in Ref. 59 

Recently, Giovannetti et al.^^ also proposed the down¬ 
folding scheme for the electron-phonon coupled systems. 
Here, we compare our scheme with that of Giovannetti et 
al. The main difference is the form of the double count¬ 
ing correction. In Ref. 59, the phonon-related part of 
the Hamiltonian,'^ which corresponds to Eq. (41) in our 
case, is given by 

'A = + ^g^P'>ni{xi - x?) 

I I 

+ (42) 

I 

Here, the electron-phonon coupling g^P'> is calculated 
at xi = x°, while in our scheme g^^^ is calculated at 
xi = 0. Furthermore, Giovannetti et al. includes the 
term '^i which represents the deformation of 

the band due to the difference in the equilibrium po¬ 
sition between that in the low-energy Hamiltonian and 
that obtained by the DFT optimization. Thus, Giovan¬ 
netti et al. also introduce the correction to the electronic 
part, while our scheme only includes the correction to 
the phonons. Therefore, the form of the Hamiltonian in 
Ref. 59 is more general than ours. 

If the difference in the equilibrium position is 
large (i.e., |x[’| is large), the band deformation term 
g^P^x^ni would become important. In the case of the 
fullerides, which will be discussed in the next section, we 
conclude that this effect is small because the equilibrium 
positions of the undoped and doped Geo solids are very 
similar, which makes the effect of doping almost rigid 
band shift. Thus, our scheme is well applicable to the 
fulleride problem. However, of course, there exist sys¬ 
tems in which this band deformation effect is significant. 
Ref. 59 argues that it is important to take into account 
the band deformation effect in the case of the K-doped 
picene system,®^ because the deformation of the molecule 
by the doping is not negligible. 

To derive the parameters in Eq. (42), Giovannetti et 
al. assume that the electrons couple to a single opti¬ 
cal phonon, while in our scheme, we can treat all the 
phonon modes. Eirst, Giovannetti et al. estimate g^^^ by 
calculating the electron-phonon coupling for the undoped 


picene system. Then, they determine and such 
that the mean-field solution of the Hamiltonian recovers 
the equilibrium positions and the phonon frequencies of 
the doped system derived within the DFT and DFPT. 

We still lack the methodology to derive the Hamil¬ 
tonian with the form of Eq. (42) in a totally ab initio 
way, i.e., without simplifying the electron-phonon cou¬ 
pling or determining and in the post processing. 
To realize this, we need to develop an ab initio struc¬ 
ture optimization scheme without the effect of the low- 
energy electrons. We also have to carefully consider the 
change of the electronic parameters. When we derive a 
model based on the optimized structure without the ef¬ 
fect of the low-energy electrons, for example, the shape of 
the Wannier function can be different from that with the 
fully optimized structure. Then, the values of Coulomb 
interaction parameters can differ from those of the con¬ 
ventional cRPA, which uses the Wannier functions con¬ 
structed from the fully optimized structure. When the 
band deformation is really severe, we might have to be 
careful in the choice of the low-energy subspace since 
the low-energy band character of the optimized struc¬ 
ture without the effect of the low-energy electrons might 
change from that of the fully optimized structure. There¬ 
fore, there left many open questions and challenges in the 
derivation of the Hamiltonian including the band defor¬ 
mation term [Eq. (42)]. Note that, in the situation where 
this band deformation is important, the cDEPT is also 
challenged, since the current cDFPT does not take ac¬ 
count of its effect. 


IV. APPLICATION 
A. Calculation conditions 

We performed the cDFPT calculations®^ for the five 
different fee A 3 C 60 systems, namely, K 3 C 60 , RbsCeo, and 
CS 3 C 60 with three different lattice parameters, whose 
properties are summarized in Table I. We employed 
the same lattice constants as those employed in Ref. 65 
to evaluate the Coulomb parameters by the cRPA. We 
specify the material by the volume occupied 

per Ceo®” anion in solid. The most expanded material 
(CS 3 C 60 with = 804 A®) is a Mott insulator and the 
second most expanded system (CS 3 C 60 with = 804 

A®) is on the verge of the metal-insulator transition.®® 
The other three materials show a metallic behavior and 
the superconductivity emerges at low temperature. 

As is already explained in Sec. HI, the implementation 
of the cDEPT can be done by slightly modifying the ex¬ 
isting DFPT program. Among the various DFPT codes, 
in the present study, we modified the one implemented 
in QUANTUM ESPRESSO package®^’®^ (see Appendix C). 
In the cDFPT calculation, we need to define the low- 
energy subspace (t-subspace). Fig. I shows the band 
structure for fee CssCeo with = 762 A®. Around 

the Fermi level, there exist the so called tiu bands origi- 


TABLE I: List of materials employed in the calculation. We 
show the name of compounds, the lattice constant a, corre¬ 
sponding volume occupied per Ceo®” anion in solid, applied 
pressure in the experiments, and superconducting transition 
temperature Tc or Neel temperature Tn. The listed materi¬ 
als are the same as those of Ref. 65 (in Ref. 65, the Coulomb 
interaction parameters and the hopping parameters were eval¬ 
uated) . 



a 

(A) 

^Ced'- 

(A®) 

Pressure 

(kbar) 

Tg (Tn) 

(K) 

Ref. 

fee K3C6O 

14.240 

722 

0 

19 

67 

fee RbsCeo 

14.420 

750 

0 

29 

67 

fee CssCeo 

14.500 

762 

7 

35 

66 

fee CssCeo 

14.640 

784 

2 

26 

66 

fee CssCeo 

14.762 

804 

0 

(2.2) 

66 



FIG. 1: (Color online) DFT band structure of fee CsaCeo with 
Vq^^- = 762 A®. We choose the tiu bands as the t-subspace 
(the yellow shaded region). The green dotted curves denote 
the band dispersion derived by Wannier hopping parameters 
for the t-subspace. 


nating from the three-fold degenerate LUMO orbitals of 
the Ceo molecule. The bands are isolated from the 
other bands. As we already mentioned in Sec. I, the 
low-energy physics is governed by the low-energy bands, 
therefore, we choose the Hu bands as target bands. 

The phonon calculations with the cDFPT and the 
DFPT were performed subsequently to the DFT ground- 
state calculations. In the DFT part, we adopted the 
local density approximation (LDA) with the Perdew- 
Zunger parameterization.®® The pseudopotentials for C, 
K, Rb, and Cs atoms were prepared with the same proce¬ 
dure as in Ref. 69 (the Troullier-Martins norm-conserving 
pseudopotentials^® in the Kleinman-Bylander represen¬ 
tation^®). We employed 4x4x4 k mesh and the cutoff 
energy of 50 Ry for the wave functions. With the above 
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conditions, we performed the structure optimization for 
the materials listed in Table I with fixing the lattice con¬ 
stant and with ignoring the orientational disorder. In the 
phonon calculation part, we employed 2 x 2 x 2 q mesh and 
the Gaussian smearing of 0.025 Ry. 


B. Phonon frequencies 

In the alkali-doped fullerides, it has been shown 
that the dominant electron-phonon coupling is coming 
from the intramolecular vibration.When we con¬ 
sider the isolated Ceo molecule, only the intramolecular 
phonon modes with the Ag and Hg symmetries have finite 
electron-phonon couplings to the tiu electrons.This 
is because the Cgo molecule has extremely high symme¬ 
try {Ifi symmetry) and the coupling to the other modes 
are forbidden due to the symmetry reason.®^ This prop¬ 
erty also holds well in the Geo solids. In particular, the 
coupling to the Jahn-Teller phonon (so called Hg modes) 
is argued to be crucial to the superconductivity.^^’®^ 

Table II summarizes our calculated partially renormal¬ 
ized phonon frequencies (w^^^’s) of the Hg modes at T 
point. Due to the crystal field, the frequencies of the Hg- 
mode are split into two. The high phonon frequencies up 
to ~ 1600cm“^ (~ 0.2 eV) can be ascribed to the stiff 
G-G bonds and the lightness of the carbon atoms. Fur¬ 
thermore, the intramolecular nature of the modes leads 
to the following features: The Hg phonon modes have 
little dispersion (see Fig. 2). The material dependence of 
the frequencies is weak. 

Note that these partially renormalized frequencies 
w^^^’s are the inputs for the low-energy solvers and thus 
can not be directly compared with the experimentally 
observed frequencies. To compare with the experiments, 
we have to include the effect of the t-subspace electrons 
and calculate the fully renormalized phonon frequencies 
(w'-'^^’s). In general, a stronger coupling between the t- 
subspace electrons and the phonons leads to a larger dif¬ 
ference between w^^^’s and In the case of the 

alkali-doped fullerides, the electron-phonon coupling of 
the individual mode is not large, while the accumulation 
of the contributions leads to the total electron-phonon 
coupling of A ~ 0.5-1.0.®“®’®^ Therefore, we do not ex¬ 
pect a large difference between w^^^’s and 

In Table III, we list the fully renormalized phonon fre¬ 
quencies of the Hg modes at F point computed by the 
DFPT. By comparing them with the partially renormal¬ 
ized values in Table II, we see the softening of the fre¬ 
quencies. This is because the phonons are dressed by the 
coupling between the phonons and the tiu electrons. In 
other words, the phonons acquire the self-energy associ¬ 
ated with the t-subspace electrons. However, as is ex¬ 
pected (see the discussion above), the difference is small: 
The absolute difference is at most ~ 50 cm“^. If we 
consider the ratio it exceeds 0.95, i.e., the dif¬ 

ference is less than 5 %. Even when we accurately treat 
the t-subspace processes beyond the DFPT level by the 
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FIG. 2: (Color online) Phonon dispersion of fee CssCeo with 
= 762 A®. For the sake of visibility, we restriet the 
frequeney range to 1100-1400 om“^. The red solid (blue dot¬ 
ted) curves indicate partially (fully) renormalized frequencies 
calculated by the cDFPT (conventional DFPT). 


model calculation method, the t-subspace renormaliza¬ 
tion effects would remain small. Then, we can expect 
that the conventional DFPT calculations give reasonable 
estimates of the phonon frequencies. Indeed, the fully- 
renormalized frequencies in Table III agree well with the 
experimental data.^^’®^ 

Figure 2 shows both the partially (red) and fully (blue) 
renormalized phonon frequencies between 1100 cm“^ and 
1400 cm“^ for fee GsaGeo with = 762 A®. Several 
intramolecular modes including Hg modes [Hg{^) and 
Hg{6)] and the non-i7g modes exist in this frequency 
range. While they are common in that they have lit¬ 
tle dispersions, we see a clear difference between the 
Hg modes and the others in the way of the softening: 
The non-Hg modes do not couple to the electrons.®® 
Hence, their frequencies are not affected by the inclusion 
of the f-subspace renormalization effects. As a result, 
the blue dotted curves (w^^^) are on top of the red solid 
curves (w^^^) for the non-Hg modes. On the other hand, 
the frequencies for the Hg modes are renormalized by a 
few percent. Indeed, the red and blue curves are located 
at different positions for the Hg modes (see the frequency 
regions 1100-1150 and 1260-1300 cm“®). 


C. Phonon-mediated effective interactions between 
the low-energy electrons 

If we write down the partition function for the electron- 
phonon coupled Hamiltonian in Eq. (15) in the coherent 
state path-integral formalism, we find that we have at 
most quadratic term for the phonon fields. Then, we 
can integrate out the phonon degrees of freedom ana¬ 
lytically. This results in an electronic model with the 
additional electron-electron interaction mediated by the 
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TABLE II: Partially renormalized phonon frequencies of Hg modes at P point calculated by cDFPT method. The unit is 
in cm“^ (1 eV = 8065.54 cm~^). In the ideal Ih symmetry (molecular limit), the TLg-mode phonon frequencies are five-fold 
degenerate. In the fee AsCeo systems, they are split into three-fold degenerate and two-fold degenerate frequencies due to the 
crystal field. Thus, we show two frequencies for each Hg mode. In the parentheses just after the material names, we show the 
values of in A®. 


mode 



frequency [cm ^ 

1 


KsCeo (722) 

RbsCeo (750) 

CsaCeo (762) 

CsaCeo (784) 

CssCeo (804) 


260, 271 

258, 269 

259, 278 

259, 274 

258, 272 


433, 435 

433, 433 

434, 436 

434, 435 

433, 435 

HA3) 

706, 708 

707, 708 

709, 710 

709, 710 

709, 710 

HgA) 

785, 786 

785, 787 

787, 797 

786, 793 

785, 791 

HA^) 

1124, 1128 

1124, 1129 

1129, 1138 

1127, 1135 

1126, 1132 

HA6) 

1282, 1287 

1282, 1287 

1292, 1298 

1288, 1294 

1286, 1291 

HA7) 

1451, 1455 

1452, 1455 

1463, 1466 

1459, 1461 

1457, 1459 

HA») 

1563, 1564 

1563, 1565 

1573, 1573 

1569, 1570 

1567, 1568 


TABLE III: Fully renormalized phonon frequencies of Hg modes at F point calculated by conventional DFPT method. The 
unit is in cm~^ (1 eV = 8065.54 cm~^). The splitting of the frequencies of each Hg mode is due to the crystal field. We show 
the values of in A® in the parentheses just after the material names. For comparison, we also show the experimentally 

observed phonon frequencies in KsCeo-^^ 


mode 



frequency 

[cm-i] 



KsCeo (722) 

RbsCgo (750) 

CsgCeo (762) CsgCgo (784) 

CssCfio (804) 

“KsCeo (expt.) 

HAA 

257, 268 

255, 267 

256, 277 

255, 273 

255, 271 

271 

HA2) 

423, 425 

422, 423 

422, 425 

421, 424 

420, 423 

431 

HAS) 

683, 686 

684, 686 

686, 688 

686, 688 

686, 687 

723 

HAA 

777, 778 

777, 778 

780, 788 

779, 785 

778, 782 


HAS) 

1110, 1114 

1110, 1114 

1116, 1125 

1113, 1121 

1112, 1118 


HAG) 

1267, 1273 

1267, 1272 

1277, 1283 

1273, 1278 

1270, 1275 


HA7) 

1402, 1407 

1403, 1405 

1415, 1415 

1410, 1410 

1406, 1407 

1408 

HAS) 

1531, 1536 

1531, 1535 

1541, 1544 

1537, 1540 

1535, 1538 

1547 


“ Raman scattering measurement. Ref. 72 



FIG. 3: Feynman diagram for phonon-mediated interaction 
between electrons. Solid lines with arrows, wavy lines, and 
dots represent the electron propagator, the phonon propaga¬ 
tor, and the electron-phonon coupling, respectively. 


phonons (Fig. 3),®"^ whose onsite (= intramolecular) part 


V^lLriiuJn) is given by 


51 








2uj, 


(p) 


qv 




~{P)* 

( P )^2 




(43) 


where is the number of q-mesh and a;„ is the bosonic 
Matsubara frequency a;„ = 2TmT with the temperature 
T.®® Here, are given by 

= (44) 

k 

Here, the partially-screened electron-phonon coupling 
is used to calculate the phonon-mediated interac¬ 
tions. In Appendix E, we discuss that the vertex correc¬ 
tion for g^P'> is small, which makes the estimate of the 
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TABLE IV: Material dependence of static part (a;n = 0) of the effective intramolecular interactions mediated by phonons. The 
values in the parentheses just after the material names denote Vcgt?- 


type of int. 



interaction [meV] 



KaCeo (722) 

RbsCeo (750) 

CssCeo (762) 

CssCeo (784) 

CssCeo (804) 


-152 

-142 

-114 

-124 

-134 


-53 

-42 

-13 

-22 

-31 

4?(0) 

-50 

-51 

-51 

-51 

-52 

U^) 

-73 

-74 

-73 

-74 

-75 


28 

29 

30 

31 

31 

4h(0) 

-51 

-52 

-52 

-52 

-53 


phonon-mediated interactions without the vertex correc¬ 
tion reliable. 

The phonon-mediated interactions Vifiiji are dy¬ 
namical interactions, which vanish in high frequency limit 
{iOn —t oo). We call the intraorbital density-density- 
type interaction, interorbital density-density-type inter¬ 
action, exchange-type interaction U'^^\iLOn), 

and respectively, i.e., 

(*Wn) = (45) 

with i ^ j. We also define the fully screened quantities 
Upl\iuJn), U'J^\iu}n), and Jpi^iuJn) in the same way, 
i.e., U^i\iu;n) = and 

Jph^(jw„) = Vl-/]^{iujn) = V^/Jj{iuJn). We find that, be¬ 
cause of the high symmetry of the tin orbitals, the val¬ 
ues of U^^^\iuJn), U'^^’^\iuJn), and do not 

depend on orbital. 

Table IV summarizes the values of the static parts of 
these interaction (a;„ = 0). We find that the relation 
^ U^h^\0)-2jl^^\0) well holds, which also 
holds for finite frequency (see Fig. 4). We first discuss 
the partially renormalized interactions. The negative val¬ 
ues of 17p((^(0), 17p|f^(0), and Jph^(O) indicate that the 
interactions are attractive at = 0. Therefore, they 
will compete with the repulsive onsite Coulomb inter¬ 
actions. As for the density-density channel, since the 
intramolecular Coulomb repulsion (the Hubbard U) for 
the tiu electrons is estimated to be on the order of ~ 1 
eV,®® the repulsive Coulomb interaction dominates over 
the phonon-mediated attraction. However, remarkably, 
the situation changes for the exchange-type interaction: 
the absolute values of |>/ph^(0)| ^ 0.05 eV is larger than 
those of the Hund’s coupling J ~ 0.035 eV.®® There¬ 
fore, in the fullerides, an effectively negative exchange 
interaction is realized.This is in constant with e.g.. 


the case of LaFeAsO (the first discovered iron-based su¬ 
perconductor®®), where the Hund’s coupling is as large 
as ^ 0.5 and the phonon-mediated exchange in¬ 

teraction Jph^O) ~ —0.02 eV gives only a minor cor¬ 
rection.®^ The unusual competition of the Hund’s cou¬ 
pling and the phonon-mediated interactions can be as¬ 
cribed, mainly, to the following two reasons.^® One is the 
molecular nature of the maximally localized Wannier or¬ 
bitals. Then, the sizes of the Wannier orbitals become 
larger than those of atomic-orbital-like Wannier func¬ 
tions, which results in a smaller Hund’s coupling. The 
other is the enhancement of the negative ( 0 ) due to 
the strong couplings between the Jahn-Teller modes and 
the tiu electrons. The Jahn-Teller Hg modes give the 
non-density-type electron-phonon coupling, which con¬ 
tribute to (0).^®’®®’®® Note that the non-Jahn-Teller 
Ag modes do not contribute, since the couplings of the 
Ag modes are of density-type. 

As for the material dependence, while that of (0) is 
small, we see discernible material dependence in 17^^ ( 0 ) 

and Hpj(^^(0). We identify the origin of the material de¬ 
pendence to be the vibration modes of the alkali ions 
at the tetrahedral sites. It is reasonable that they give 
a material-dependent contribution as the distances be¬ 
tween the Ceo” anions and/or the alkali cations change. 
Indeed, if we compute (0) and C/p|f^(0) for the five 
materials with excluding the alkali-ion contributions (in 
this case, the values become the sum of the contribution 
from the intramolecular phonons), the results for ( 0 ) 
[C/pj(’^(0)] are, in ascending order of , —89 [10], —91 
[9], —91 [9], —93 [ 8 ], and —95 [ 8 ] meV. As is clear, they 
have much less material dependence than those with the 
alkali-ion contributions, which is natural because we can 
expect that the intramolecular phonons have little mate¬ 
rial dependence. We find that the alkali-ion modes cou¬ 
ple to the total density of the electrons, i.e., they 
couple to the density of the individual orbital with al¬ 
most the same amplitudes ~ ~ ^ 33 ^ )■ Thus, it 
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(b) imaginary part, real frequency 
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FIG. 4: (Color online) Frequency dependence of phonon-mediated interactions for fee CssCeo with = 762 A® along 

the real frequency [panels (a) and (b)], and along the Matsubara axis [panel (c)[. (a) [(b)] Real [imaginary] part of U^^{u!), 
and The frequency dependence of these quantities are calculated at uj + iri with rj = 0.01 eV. (c) Real part of 

and The imaginary part is always zero along the Matsubara 

axis. 


does not contribute to . The contribution to Jp^ (0) 
originates from the intramolecular Jahn-Teller coupling 
(coupling to the Hg modes). Therefore, we see little ma¬ 
terial dependence in Tp^^(O). Since the electron-phonon 
coupling of the alkali-ion modes are of density-type, the 
alkali-ion mode contribution is efficiently screened by the 
flu electrons, which leads to a minor role of the alkali- 
ion modes in the superconductivity. As a result, as we 
will see below, the dominant contribution to the fully 
renormalized interactions comes from the intramolecu¬ 
lar phonons, which is consistent with the previous stud¬ 
ies. 

We can compute the fully-screened phonon-mediated 
onsite interactions using Eq. (43) by replacing the par¬ 
tially renormalized quantities with the fully renormal¬ 
ized quantities. We list the values of their static part 
(ujn = 0) in Table IV. We find that the magnitudes of 
density-density type interactions, C7p/^(0) and C7p^^(0), 
differ substantially from those of the partially renormal¬ 
ized ones, and 17p^^(0). On the other hand, the 

values of TpP(O) are almost unchanged from those of 

Jp^^(O). This different behavior between 17ph, C7ph and 
Jph can be understood as follows. The tiu electrons effi¬ 
ciently screen the non-Jahn-Teller type electron-phonon 
coupling, while the Jahn-Teller type coupling not. The 
former contributes to Uph and Therefore, the differ¬ 


ence between the partially and fully renormalized quan¬ 
tities is substantial. On the other hand, only the Jahn- 
Teller phonon contributes to Jph. Therefore, we have lit¬ 
tle difference between the partially and fully renormalized 
quantities. As we discuss above, the alkali-ion mode con¬ 
tribution becomes small in the fully renormalized quan¬ 
tities and the intramolecular Hg mode contribution be¬ 
comes dominant (intramolecular Ag mode contribution 
is also screened because Ag mode couples to the total 
density of the electrons), which makes the material 
dependence of Up(^{0), Up^\o), and Jph^O) small. 

When we consider the contribution from Hg modes in 
the molecular limit,we can show that the relation 
Up-l\iuJn) = —Up(^{iLUn)/2 holds. Since, in reality, we 
have a small contribution from the other modes such as 
Ag and the alkali-ion modes, the above relation does not 
exactly hold. However, this naturally explains why the 
interorbital interactions become repulsive (Up/^O) > 0). 

We finally discuss the frequency dependence of the 
phonon-mediated interactions. The frequency depen¬ 
dences for fee CsaCeo with = 762 on the real 

frequency axis are shown in Figs. 4(a) and (b), where 
the panel (a) [(b)] shows the real [imaginary] part of the 
phonon-mediated interactions. Since the frequencies of 
the intramolecular phonons lies up to ^ 0.2 eV, there ex¬ 
ist significant structures below ^ 0.2 eV. Im U^ioj) and 
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Im (w) are always negative. On the other hand, the 

Im can be both negative and positive. This is 

because the contribution from the non-Jahn-Teller and 
the Jahn-Teller {Hg) phonons coexist [the former (lat¬ 
ter) gives negative (positive) contribution]. Note that 
both the non-Jahn-Teller and Jahn-Teller phonons give 
negative contribution to Im (uj) and Im 

We also show the frequency dependence along the 
Matsubara frequencies in Fig. 4(c), where we also plot 
the frequency dependence of the fully screened interac¬ 
tions. Again, the non-monotonic behavior in U^^\iu}n) 
can be ascribed to the coexistence of the attractive 
(non-Jahn-Teller) and repulsive (Jahn-Teller) contribu¬ 
tions. Because the contribution from the non-Jahn-Teller 
phonons becomes small in the fully screened interactions, 
the frequency dependence of Up\^\iuJn) becomes mono¬ 
tonic. Since only the Jahn-Teller modes, which is poorly 
screened by the tiu electrons, contribute to Jph, we have 
small difference between (iuJn) and Jp{^(iaj„). Fi¬ 
nally, we note that the relation 

well holds along both the real and imaginary frequency 
axes. 


In this paper, we have focused on the alkali-doped ful¬ 
lerides. However, in principle, the cDFPT is applicable 
to other materials in which phonons play a crucial role. 
These applications remain as interesting and important 
future issues. There also remain challenges in the de¬ 
velopment side: As discussed in Sec. HIE, the current 
cDFPT is not applicable when the equilibrium positions 
of the ions change drastically by the coupling to the low- 
energy electrons. It is also challenged when there exists 
strong anharmonicity in the system. These are important 
open questions in the downfolding for electron-phonon 
coupled systems. 
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V. CONCLUSION AND OUTLOOK 


Appendix A: Equivalence of Eqs. (6) and (8) 


In this paper, we have presented a detailed explana¬ 
tion on the newly developed ab initio downfolding scheme 
for the electron-phonon coupled system, cDFPT. With 
the cDFPT, we can calculate the partially renormalized 
phonon frequencies and electron-phonon coupling, which 
is used as the parameters in the effective low-energy 
Hamiltonian. We have shown that the cDFPT scheme 
can be easily implemented by the slight modification of 
the conventional DFPT scheme. 

We have applied the cDFPT scheme to the alkali- 
doped fullerides. By excluding the t-subspace renormal¬ 
ization effect, we have seen the hardening of the frequen¬ 
cies of the phonon modes which couple to the t-subspace 
electrons. We have also discussed the difference be¬ 
tween the partially and fully screened phonon-mediated 
interactions. In the partially screened interactions, the 
non-Jahn-Teller phonons give substantial contributions. 
However, in the fully screened screened interactions, the 
contribution from the non-Jahn-Teller modes becomes 
small because it is efficiently screened by the elec¬ 
trons. Then, the Jahn-Teller phonons give the dominant 
contributions to the fully screened interactions. 


Here, we show that Eqs. (6) and (8) indeed give the 
same solution. When we write Eq. (8) as 

(^SCF + Q — £n) |A^„) = —{Of.ti — T’n)AI/sCF|'0n), 

AX y 

(Al) 


the A matrix is given, in the Bloch basis (note that, in 
the case of quantum espresso, the plane basis is used 
in the actual calculation), by 


A = 


£l -|-Q!i — £„ 


0 \ 


£ 2 -\- 0-2 — 


0 


£M-\-C»-M — £n j 


(A2) 


where M is the size of the basis set to describe the Bloch 
states. — Pn on the r.h.s of Eq. (8) is rewritten as 
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OF,n -Pn=J2 

m 

= E 


^F,n(l ^n,77i) ^F,m^m,n ' 


^F,n — ^F,r 

i 


|‘0m){‘07i 


/a a \ a ^F,m 2 

(^L'F.n ^F,m)^m,n (^m ^m,n 




E ^F,n — ^F,m ^ \ 

“r 0!.jji Sjij 

Sn £--- 


n ^m 


V^m) {V^m I • 


(A3) 


With Eqs. (Al), (A2), and (A3), we can show that | Ai/)n) 
is given by 

lA-t/jn) = A~^y 

= ^ _ ~~^’"" ^m,rt|t/’m)(V^m| AVsCfI'^ti), 

(A4) 

which is nothing but a proof that Eq. (8) gives the same 
result as that of Eq. (6). 

Appendix B: Comparison between cDFPT and 
cRPA 

Here, we compare the present cDFPT with the 
cRPA.^^ In the cRPA, which derives the effective 
electron-electron interactions in the low-energy model, 
we calculate the partially screened Coulomb interaction 
as23 

= { l - VXr )~^ V . (Bl) 

The fully screened Coulomb interaction is obtained by 
further taking into account the t-subspace screening ef¬ 
fect: 

wif) = (l-W‘'P'^xf)~^ W'Py (B2) 

One can see that Eqs. (Bl) and (B2) have the same 
structure as that of the screened electron-phonon cou¬ 
pling [Eqs. (30) and (31)]. Both the cRPA and cDFPT 


methods relies on the same kind of decomposition of the 
screening processes. In both cases, we calculate the par¬ 
tially screened quantities, which is to be used in the low- 
energy Hamiltonian. 


Appendix C: Practical implementation in the case of 

QUANTUM ESPRESSO 


Here, we provide an example how we modify a source 
code. In the DEPT implemented in the version 4.3.1 of 
QUANTUM ESPRESSO, the parameters are de¬ 

fined in “orthogonalize.fOO”, which exists in “PH” folder. 
In Ref. 90, we distribute a modified “orthogonalize.f90” 
under the GNU General Public License.®^ 

Appendix D: Confirmation of the equality 
E = Et -I- Er in Sec. Ill B 


Here, we show that the equality E = E( -|- in 
Sec. HIB indeed holds. In principle, the self-energy E, 
the electron-phonon coupling g, the polarization function 
X®, and so on, are expressed as matrices. In this section, 
for the sake of simplicity, we treat them as if they were 
scalar quantities. One can easily extend the proof to the 
case where they are matrices. Et = Is^^^Pxdft i® rewrit¬ 
ten as 
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Xt 


1 - W(p^Xt 

= \qib)\2 - 1 - 

1 - {ix° 1 - Vl^(P)x? 1 - «Xr 


Xt 


VXr 


= (^1 + Y 


Xt 


-vx°J 1-T^(p)x° 


1 + 


VXr 


Xt 


Xt 


1 - W^p'lxt 

n W^P'^ 


1 - vx° 

(i + iy(p)x°) 


■Xr 


x? + x?- 


^(p) 




1 - W(p)xt 1 - W^^P^X? 1 - W(p')x't 

Xt + XtW^^ht + XrW^^ht + XtW^^hr + 


oX° + X°W^^^^ 


Xt 


1 - W^P^Xt 




Xt 


1 - w(p)x°t 




(Dl) 


Similarly, Xdpx is rewritten as 




2 Xr 


^-VXt 
Using the equality 






Xr+Xr 


(D2) 


w(p) + w^p) _ ^ _ wip) =__= wU) 

1 - iu(p)x? 1 - W(p)x°t 


(D3) 


one can show that Et + is expressed as 
Et + E,=|g('') \^\x°t+Xr + {Xt + Xr) {Xt + Xr) 




X° + 

= |5^^^Pxdft, (D4) 

which agrees with the expression for E in Eq. (33). 


Appendix E: Smallness of electron-phonon vertex 
correction in downfolding procedure 

Due to the high phonon frequency ~ 0.1 eV, which is 
comparable to the typical electronic kinetic energy ^0.5 
eV, the Migdal theorem®^ is violated in the A 3 Geo sys¬ 
tems. Therefore, we need a careful consideration about 
the vertex corrections. In this Appendix, we argue that, 
as far as the processes involving the high-energy electrons 
are concerned, the electron-phonon vertex corrections are 
small. 

To see this, let us consider the “first-order” vertex cor¬ 
rection diagram in Fig. 5. For simplicity, we assume 
that the multiple intramolecular phonon modes are rep¬ 
resented by a single Einstein phonon branch with the 


frequency ujq and that the electron-phonon vertex g has 
no momentum dependence (or the electron-phonon cou¬ 
pling is local). Then, the inclusion of the diagram in 
Fig. 5 gives the correction to the bare electron-phonon 
vertex go as go —>■ 50(1 +7) with 7 being a dimensionless 
quantity given by 

7 = ^ gig 2 D{k - k')G{k')G{k'+ g), (El) 

k' 

where g^’s, D, and G are the dressed electron-phonon 
vertices, phonon Green’s function and electron Green’s 
function respectively, gt’s and G have orbital indices, 
while we do not show them for simplicity, k [q] repre¬ 
sents a set of the momentum and the fermionic [bosonic] 
Matubara frequency k = (k,z/ji) [q = (q, a;„)]. T is the 
temperature and Ak is the number of k-points. Note 
that this diagram is of first order with respect to D, 
however, it contains the higher order diagrams with re¬ 
spect to the bare phonon Green’s function Dq- In the 
downfolding procedure, the low-energy processes are ex¬ 
cluded, therefore, the two electron Green’s function in 
Eq. (El) should be a combination of Gh and Gh or of 
Gh and Gl, where Gh {Gl) is the propagator of the 
high-energy (low-energy) electrons®^. Then, the typical 
order of 7 associated with the downfolding is given by 
I 7 I ~ 2gig2/u}r X l/AE with the renormalized phonon 
frequency ujr and the typical particle-hole excitation en¬ 
ergy scale involving high-energy degrees of freedom AE. 
Here, to derive this expression, we have employed the fact 
that the typical order of the convolution of Gh and Gh 
or Gh and Gl is ^ l/AE. In the case of the alkali- 
doped fullerides, AE is at least ^ 1 eV. 2 ^ 132 /wr is 
nothing but the static part of the fully-screened phonon- 
mediated interaction. If gi and g 2 are the coupling be¬ 
tween phonons and low-energy electrons, it corresponds 
to Uli\0), U'J^\0) and J^P(O) in Table. IV. While we 
do not estimate the coupling between the phonon and the 
high-energy electrons, we expect the order of the phonon- 
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FIG. 5: First-order electron-phonon vertex correction dia¬ 
gram. The open (closed) circle represents the bare (dressed) 
electron-phonon coupling. The bold (thin) solid and wavy 
lines indicate the dressed (bare) electron and phonon propa¬ 
gators, respectively. 


mediated interactions involving the high-energy electrons 
is the same as that of [/pjf^(O) and Jp{^(0). In addition, 
in the diagrams considered in the downfolding procedure, 
the orbital indices for gi and g 2 are usually different. This 
is because one is the coupling to electron state and the 
other is the coupling to the hole state, which would make 
2gig2loJr smaller. In any case, 2g\g2l^r will be at most 
^0.1 eV. As a result, the correction 7 associated with 
the downfolding will take a small value 7 < 0 . 1 . 

In conclusion, the neglect of the electron-phonon ver¬ 
tex corrections in the model-derivation step as in the case 
of the cDFPT is justified. However, we note that we still 
need a careful treatment for the vertex corrections in the 
model-analysis step since the vertex corrections in the 
t-subspace is not negligible any more. 
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